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ABSTRACT 



Context. Over the past decade the use of the red clump (RC) as distance indicator has increased in importance as this evolutionary 
phase is well populated and a good local calibration exists. 

Aims. The absolute calibration of the RC in the / and K band is investigated again based on the recently published revised Hipparcos 
parallaxes. 

Methods. A numerical model is developed that takes the various selection criteria and the properties of the Hipparcos catalogue into 
account. The biases involved in applying certain selections are estimated with this model. 

Results. The absolute magnitudes that are found are Mi = -0.22 ± 0.03 and Mr = -1.54 ± 0.04 (on the 2MASS system). The 
7-band value is in good agreement with previous determinations, the Tf-band value is fainter than previously quoted, and this seems to 
be related to a selection bias whereby accurate Tf-magnitudes are only available for relatively few bright stars. Applying population 
corrections to the absolute K magnitude of RC stars in clusters supports the fainter magnitude scale. 

Key words. Stars: distances - distance scale 
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1. Introduction 

Paczyriski & Stanek (1998) were the first to advocate the use of 
red clump (RC) stars in the /-band as distance indicators. Two 
advantages are that this phase of stellar evolution is usually well 
populated, and that many RC stars exist with good parallaxes in 
the Hipparcos catalogue. 

They consider stars with an error in parallax smaller than 
10%, located in the box 0.8 < {V-I) < 1.25 and -1.4 < Mi < 
1.2 and with Hipparcos flag H42 equal to one of "A,C,E,F,G", 
which means that the star has one or more /-band measure- 
ments in the Cousins, Johnson, or Kron-Eggen system. They 
find no significant trend of Mi on colour and derived Mi — 
—0.185 ± 0.016 for the entire sample. They then considered the 
228 stars with d <70 pc and average distance of 50 pc to find 
Ml = -0.192 ± 0.023, and A/j = -0.094 ± 0.027 for the rest 
of the sample at an average distance of 106 pc. The difference 
is ascribed to reddening, and when extrapolating to pc, they 
arrive at their final result of find Mi = —0.279 and assign an 
error bar of 0.088. In deriving these number, they also applied a 
correction for "a distance bias". 

Stanek & Garnavich (1998) considered the same sample, and 
their final value of Mi = —0.23 ± 0.03 is based on the sub 
sample of stars within 70 pc and they argue that no coiTection 
for reddening is needed, and they do not apply a correction for 
selecting the sample on parallax. 

Udalski (2000) considered the effect of metallicity. He took 
the same cut in parallax error and Hipparcos H42 flags, but 
also included stars with flag "H" to increase the sample and 
considered stars with a metallicity determination from high- 
resolution spectroscopy by McWillam (1990). Based on this 
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sample of 284 stars, he finds a weak dependence on metallic- 
ity: Mi = (0.13 ± 0.07)([Fe/H]+0.25) + (-0.26 ± 0.02). 

Alves (2000) was the first to consider the /C-band for RC 
stars as a standard candle. He considered the 284 stars from 
Udalski (2000) and finds that 238 of them have a iiT-magnitude 
in the Two Micron Sky Survey (TMSS, Neugebauer & Leighton 
1968; stars are designated by "IRC" numbers). Assuming no red- 
dening and considering stars with —2.5 < Mk < —0.8 and all 
{V — K) colours, he finds Mk = — 1.61 ± 0.03 and a weak de- 
pendence on metallicity Mk = (0.57±0.36)[Fe/H]+(-1.64± 
0.07). 

The RC has gained importance as a distance indicator since 
the Paczyriski & Stanek (1998) paper as shown for example by 
distance estimates based on the RC to Fornax (Rizzi et al. 2007), 
M33 (Kim et al. 2002), or IC 1613 (Dolphin et al. 2001). 

With the recent publication of the revised Hipparcos parallax 
data, it is timely to investigate the absolute magnitude of the 
RC again. In Sect. 2 the revised Hipparcos data is introduced 
with auxiliary data used in the analysis. Section 3 describes the 
numerical model used for estimating selection biases. Section 4 
describes the results, and Section 5 presents a brief summary and 
discussion. 



2. The Data 

The Hipparcos data used are the parallax and error in the parallax 
from the re-reduction of the raw data by van Leeuwen (2007). 
For other data {V, /, Hipparcos flags), the originally published 
catalogue is used (ESA 1997). 

Using the script submission tool within VizieR (Ochsenbein 
et al. 2000) the list of Hipparcos objects was correlated with 
selected measurements available, in particular, spectroscopic 
[Fe/H] values (available for 2445 stars, and adding 11 stars 
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from Zhao et al. 2001, 64 from Liu et al. 2001, and 178 from 
Mishenina et al. 2006), the TMSS catalogue (Neugebauer & 
Leighton 1969, available for 3167 stars), and Stromgren pho- 
tometry (available for 36337 stars). In case of multiple [Fe/H] 
determinations the mean was taken. 

Independently, the Hipparcos objects are correlated on posi- 
tion with the DENIS (Epchtein et al. 1999, for 42696 stars) and 
2MASS (Skrutskie et al. 2006, for 1 17358 stars) NIR catalogues. 
The TMSS was the main source of i^T-magnitudes used by Alves 
(2000). This is the only source of uniform /^-magnitudes for 
bright objects that often saturate in the DENIS and 2MASS sur- 
veys. In the present study TMSS, DENIS and 2MASS are be- 
ing used but this requires transformation to the 2MASS system, 
which will be the system of choice. Using 749 objects with a 
2MASS photometric quality flag, phqual-flag, of "A" in the K- 
band and that not saturate in DENIS / (implying / > 10) it 
is found that DENIS K - 2MASS K = -0.018 ± 0.066 with 
no colour term on (/ — i^)DENis- There are no stars in com- 
mon between the TMSS survey and non-saturated DENIS or 
2MASS objects with phqual-flags of "A" or "B" in K. The trans- 
formation to the 2MASS system is done indirectly via the SAAO 
Carter (1990) system. The advantage over a comparison to the 
Koornneef (1983) system, which was used by Alves (2000), is 
that the colour term and offset with respect to the 2MASS system 
are smaller (Carpenter 2001, and updated on the IPAC websit^H). 

Using 62 TMSS stars that are listed in Carter (1990) and that 
have both an IRC K and I magnitude available, the following 
relations are derived: (J — iir)cartcr — 0.484(7 — K)mc ^ 
0.103 with a dispersion of 0.04 mag, and eliminating 5 outliers, 
-fi^Caitcr — Kijic — —0.011 with no significant colour term, 
and with a dispersion of 0.04 mag. The transformation from the 
Carter to the 2MASS system follows the relation on the quoted 

IPAC website: Ksmass = i^Cartcr -0.024 + 0.017( J-7C)cartcr. 

The adopted if-band magnitude is, in order of preference, 
(1) the 2MASS value, but only when phqual is "A", (2) The 
DENIS value (for non-saturated stars with /dcnis > 10) cor- 
rected by H-0.018 mag, (3) the TMASS/IRC value, transformed 
to the 2MASS system as described above. If no /jrc is available 
a mean colour for RC stars of (/ — K)m^c = 1-55 is used in the 
transformation. 

Of all Hipparcos stars that have a /T-magnitude assigned, 
3116 are based on TMSS, 17 on Denis and 103827 on 2MASS. 



3. The model 

In this section a model is described to construct synthetic sam- 
ples of stars and apply various selection criteria in order to gener- 
ate samples of stars that closely resemble in nature the observed 
data. The model is follows largely Groenewegen & Oudmaijer 
(2000) but is based on the properties of the revised Hipparcos 
catalogue. 

3.1. Galactic model 

The coordinate system used is cylindrical coordinates centred on 
the Galactic centre. The galactic distribution is assumed to be a 
double exponential disk with a scale height H in the z-direction 
(the coordinate perpendicular to the galactic plane), and a scale 
length Rqc in the galacto-centric direction. The distance to the 
Galactic centre is taken to be 7800 pc (Zucker et al. 2006). 



' www.ipac.caltech.edu/2mass/releases/allsky/doc/sec6_4b.html 



Table 1. Comparison of different reddening models 





Drimmel et al. 


Parenago 


Arenou et al. 


Marshall et al. 


1.01 ±0.82 


1.00 ±0.88 


0.99 ± 0.76 


Arenou et al. 


1.20 ± 0.62 


1.15 ±0.55 




Parenago 


1.06 ±0.24 







3.2. Luminosity function 

The luminosity function can be arbitrary but, in the case of sim- 
ulating RC stars, is assumed to be a Gaussian. In / with a mean 
Ml = —0.27 and a of 0.20. The {V — I) colours is assumed to 
be a Gaussian with mean 0.98 and <t = 0.085. In K the Gaussian 
has a mean Mk = —1.60 and a of 0.22, while the {V — K) 
colours is taken as a Gaussian with mean 2.32 and a of 0.21. 

3.3. Reddening 

Visual extinction in the simulation and to de-redden the observa- 
tions is based on several 3-dimensional reddening models avail- 
able in the literature. Marshall et al. (2006) presents reddening 
in the i^-band with 15' sampling along 64 000 lines-of-sight in 
the direction | / |< lOOdegand | 6 |< lOdeg based on 2MASS 
data. Drimmel et al. (2003) presents a reddening model based on 
the dust distribution model of Drimmel & Spergel (2001), which 
is based on COBE/DIRBE data. Arenou et al. (1992) presents a 
reddening model based on a comparison of observed B, V pho- 
tometry with predicted photometry from a star's spectral type. 
Finally, a simple law following Parenago (1940) was considered: 

Aw = 0.097/ I sin 6 I (1.0 - exp(-0.0111 d | sin 6 |)) (1) 

The Marshall et al. results have been taken as reference model, 
and their reddenings were transformed to Ay as ylK/0.12. 

However, this model is only available over a limited range 
in galactic coordinates, therefore a Monte Carlo simulation was 
made to randomly generate galactic coordinates and distances 
assuming a constant number density of stars within 3 kpc dis- 
tance of the Sun. If the galactic coordinates were within the 
range of applicability of the Marshall et al. model, the visual 
extinction from the other models was determined. In the end, 
using several 1000 positions, the average and dispersion in the 
various ratios between extinction models were calculated, and 
the results are in Table [T] in the sense model listed in the row 
divided by model listed in the column. 

Although the dispersion in all the ratios is quite large, they 
mostly agree in the mean. It only seems that the Arenou et 
al. model gives higher extinction than the Drimmel et al. and 
Parenago model. The finally adopted visual extinction was the 
average of the models by Parenago, Drimmel et al., and Arenou 
et al. multiplied by 0.84. 

3.4. Properties of Hipparcos data 

The adopted completeness function in Hp magnitude is 

(1.0 + cxp((i/p - 8.695)/0.5206))"^-^^'* . (2) 

This was derived by comparing the observed magnitude distri- 
bution to that predicted with the TRJLEGAL Galactic model 
(Girardi et al. 2005). 

The completeness functions quoted below are derived by in- 
specting the ratio of stars that fulfil a particular selection to all 
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stars, as a function of Hp magnitude. The completeness function 
of stars with /-band flags "A+C+E+F+G" is approximated as 

0.90 Hp < 5 

0.50 - 0.45 (Hp - 5)/2 5 < Hp < 7 (3) 
0.05 Hp > 7. 

If stars with /-band flag "H" are also included, this becomes 



0.92 

0.60- 

0.20 



0.40 (Hp - 5)/4 



Hp < 5 
5 < Hp < 9 
Hp > 9. 



(4) 



As outlined above. K-hand magnitudes were to Hipparcos 
stars based on high-quality 2MASS and IRC data. The complete- 
ness function is approximated as 



0.60 
0.69 

0.69 + 0.30 (Hp ■ 
0.99 



7)/1.5 



Hp < 5 
5 < Hp < 7 
7 < Hp < 8.5 
Hp > 8.5. 



(5) 



For fainter objects, a (reliable) /^-magnitude is available for al- 
most all objects because of the all-sky nature of the 2MASS sur- 
vey, but for brighter stars this is clearly not the case. 

As outlined above, spectroscopic [Fe/H] values have been 
retrieved using Vizier and recent literature. The completeness 
function is approximated as 



0.55 
0.55- 
0.02 
0.0 



0.53 (Hp - 4)/3 



Hp<4 
4 < Hp < 7 
Hp > 7 
Hp > 8. 



(6) 



The error on the Hp magnitude is small but has nevertheless 
been taken into account. The error is Gaussian distributed with 
the sigma value 

log cTHp = 0.290Hp - 5.510 Hp > 11 

logo-Hp = 0.177//J9- 4.265 Hp < 11. (7) 

As explained in the explanatory supplement to the Hipparcos 
catalogue, stars with an ecliptic latitude more than 47 degree 
were observed more often, and this resulted in lower parallax 
errors on average. When /3 > 47° the error on the parallax (in 
mas) is in the mean 

logcr^ = 0.1610//P- 1.540, 

with dispersion 0.075, and else 

log cr^ = G.1573//P- 1.368, (8) 
with dispersion 0.063. The minimum error on the parallax is 

(Iogf7,)min = 0.1610//J5-1.7 (9) 

In the simulation a parallax error is generated using the mean re- 
lation and the Gaussian dispersion. If this parallax error is above 
the minimum allowed value the star is retained. 

The "measured" parallax is based on the true parallax (from 
the true distance) and the Gaussian distributed parallax error. 



3.5. The model in practice 

To begin with, a number of input parameters need to be set: 

- scale length and scale height of the population, 

- mean magnitude and dispersion in the luminosity function in 
the reference magnitude, A (i.e. I or K in the present paper), 

- mean magnitude and dispersion in the luminosity function in 
the {V — A)-colour, 

- A\/Ay value, and 

- the number of simulations, and number of stars per simula- 
tion. 

The consecutive steps in the model are the following: 

- Three random numbers are drawn to select the distance to 
the Galactic plane, the distance to the Galactic centre and 
a random angle (j> between and 2tt in the Galactic plane 
centred on the Galactic centre. From this the distance d to 
the Sun is calculated. The galactic coordinates and ecliptic 
latitude is determined. 

From the properties of the observed sample (e.g. from the ob- 
served parallax, or a photometrically determined parallax), it 
is already possible to eliminate stars beyond a certain limit- 
ing distance. 

- The reddening in the V-band and in the reference magnitude 
is determined. 

- The true absolute magnitude in V and in the reference mag- 
nitude is determined. 

- The true Hp magnitude is determined, from the distance, 
Mv, Ay, and the fact that {Hp - V) is typicafly H-0.1 (vol. 
1, sect 1.3, p. 59 of Hipparcos explanatory catalogue). 

- The error in the Hp magnitude is determined and the ob- 
served Hp determined. 

- The completeness function in Hp is determined, and a ran- 
dom number is generated to evaluate if the star is "observed". 

- Depending on the selection of the sample, other complete- 
ness functions are evaluated and applied. 

- If the star is still "in", Eqs. (8-9) are evaluated to determine 
the parallax errors, and the observed parallax. 

- Depending on the selection of the sample, additional selec- 
tion criteria on parallax or (relative) parallax error are ap- 
plied. 

- The observed magnitude in the reference colour is deter- 
mined. 

- The output of the model typically is the "observed" paral- 
lax, error in the parallax, true distance, photometric paral- 
lax based on the adopted absolute mean magnitude. Hp, ob- 
served reference magnitude. Ay and galactic coordinates. 

4. The red clump 

Table|2]lists the observed absolute magnitudes depending on var- 
ious selections. These selections largely follow the selections 
made by previous works as quoted in the Introduction. In all 
cases, only stars where Hipparcos flag H30 (the goodness-of-fit) 
is <5, TT > 0, and single stars (isoln = 5 as listed in the revised 
Hipparcos catalogue) are considered. 

Following earlier work, stars are selected inside the box (all 
colours and magnitudes are understood to be de-reddened ac- 
cording to the model described in Sect. 13.31 ) 0.75 < — /) < 
1.25, —1.6 < A/i < 1.2. The data is binned in absolute magni- 
tude and then fitted with the formula introduced by Paczynski & 
Stanek(1998): 

N = ai + a2{m — mg) + 03(771 — mo)^ 
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exp(--((TO-mo)/aRc)^). (10) 



^RC 



The results in the A'-band are based on a subset of the se- 
lected stars in the /-band. In this case Eq. [TO] is fitted to the 
stars inside the box defined by 1.8 < iy ~ K) <2.8 and 
—3.1 < Mk < —0.1. The location of the boxes in a colour- 
magnitude diagram are illustrated in Fig. [T] An illustration of 
the fit of Eq. [TO]to the data is shown in Fig. |2] 

Several of the entries in Table|2]are put between parentheses, 
to indicate a lower quality. In the /-band this refers to selections 
that include lower-quality {y ~T) data. In K this refers to model 
5 that includes too few stars for an accurate fit of the Gaussian. 

As the selection box and Eq [TO] are designed to include a 
"background" of stars, the fit of the absolute magnitude against 
colour or metallicity is made for stars that are located inside a 
smaller box. Based on the fit with Eq. [TO] and the coefficients 
04, mo, ctrc- the number of red clump stars can be estimated. 
The location of the smaller box is chosen such that it approx- 
imately includes this number of stars. In the /-band this is 0.8 
< (V^ - /) < 1.2 and -0.6 < A/i < 0.1. In the /sT-band this is 
1.9 <{y ~K)< in and -1.9 < Mk < -1.35. 

First the dependence on colour and metallicity is investi- 
gated, fitting a linear relation. When Mi is fitted against {y — /) 
colour, the slope is 0.27 ± 0.11, 0.10 ± 0.19, and -0.02 ± 0.22, 
respectively, when all, "ACEFGH" and "ACEFG" Hipparcos /- 
band flags are considered. For these three selections, the slope 
is 0.08 ± 0.06, 0.06 ± 0.07, and 0.09 ± 0.08 when A/j is fit- 
ted against [Fe/H]. The three determinations are consistent, as 
one would expect since any dependence on metallicity should 
not depend on the origin of the /-band data. No significant de- 
pendence on colour is found, in accordance with Paczynski & 
Stanek (1998). Also the dependence on metallicity is marginal. 

When A/k is fitted against {y — K) colour, the slope is 
-0.15 ± 0.07, -0.15 ± 0.08, and -0.13 ± 0.09, respectively, 
when all, "ACEFGH" and "ACEFG" Hipparcos /-band flags are 
considered. Fitted against [Fe/H], the slopes become, respec- 
tively, 0.02 ± 0.07, -0.03 ± 0.08, and -0.04 ± 0.09. In K there 
is therefore no evidence of any metallicity dependence and only 
a marginal one on {V — K) colour. 

For the selections in Table |2l the results of the simulations 
are listed in Table [3] As mentioned in Sect. 3.2, the input values 
are A/j = -0.27 with a = 0.20, and Mk = -1.60 with a = 
0.22. One can observe that the bias can be up to 0.07 mag in 
/ depending on the selection bias. The fifth column in Table [3] 
lists the "true" value, which is the observed value from Tabled 
corrected by the bias, estimated from the difference between the 
input value and the output value from the simulation. 

The average of the 4 available models or the one with small- 
est correction (model 5) gives the identical result: Ml = — 0.22± 
0.03. In K the biases are larger, up to 0.1 mag, and the average 
of 7 estimates is A/k = -1.54 ± 0.04. 



5. Summary and discussion 

Using revised Hipparcos parallaxes and extensive numerical 
simulations that take into account de properties of the Hipparcos 
catalogue in terms of completeness and parallax error, and also 
taking into account the selection functions that are involved in 
selecting stars with certain properties (/f-band, high-quality /- 
band, spectroscopic metallicities), the absolute magnitude of the 
red clump in / and K is investigated. 

In the /-band, the value of Mi = -0.22 ± 0.03 essentially 
agrees with the value quoted by Stanek & Garnavich (1998). 



a o 



0.5 



1.5 



(V-I)o 



1.5 



2.5 



(V-K)o 



Fig. 1. Colour-magnitude diagrams for model 3. The outer box 
is used to select the data to fit Eq.[TO]to. The inner box is used 
to fit the dependence of the absolute magnitude against colour or 
metallicity. 



There is no dependence of Mi on {V — /) colour, and the depen- 
dence on metaUicity is marginal: (0.08 ± 0.07) ([Fe/H] + 0.15). 

In the /C-band, the influence of selection biases appears 
stronger and the value derived A/k = —1.54 ± 0.04 differs 
significantly from the value originally derived by Alves (2000), 
which is almost 0. 1 mag brighter, when his value of — 1 .61 on the 
Bessefl & Brett (1988) system is converted to the 2MASS sys- 
tem (Carpenter 2001). There is no dependence of A/k on metal- 
licity, and the dependence on colour is weak: (—0.15 ± 0.07) 
{{V - K)o - 2.32). 

The RC stars in clusters have also been used to derive A/k, 
originally by Grocholski & Sarajedini (2002). Absolute magni- 
tudes are derived by combining the observed mean magnitude 
in a box in colour and magnitude with a reddening estimate 
and main-sequence fitting distances. Grocholski & Sarajedini 
(2002) used the second incremental data release of 2MASS and 
14 clusters to find an absolute /^-magnitude that is in agree- 
ment with the value in Alves (2000). Using the All-Sky data re- 
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odel 






/-flags 


[Fe/H] 


N 


Ml 


CTRC,! 


N 


Mk 


0"RC,K 


1 


< 0.1 




all 


No 


3395 


(-0.206 ±0.011) 


0.249 ±0.013 


1593 


-1.398 ±0.017 


0.369 ± 0.021 


2 


< 0.1 




ACEFGH 


No 


1411 


(-0.245 ± 0.014) 


0.217 ±0.016 


436 


-1.534 ±0.029 


0.277 ± 0.036 


3 


< 0.1 




ACEFG 


No 


795 


-0.269 ± 0.016 


0.205 ±0.018 


236 


-1.582 ±0.035 


0.271 ± 0.043 


4 


< 0.1, TT 


>14.3 


ACEFG 


No 


182 


-0.259 ± 0.023 


0.193 ±0.026 


116 


-1.624 ±0.025 


0.146 ±0.025 


5 


< 0.1, TT 


>20.0 


ACEFG 


No 


86 


-0.232 ± 0.033 


0.184 ±0.035 


60 


(-1.628 ±0.033) 


0.091 ± 0.029 


6 


< 0.1 




all 


Yes 


487 


(-0.217 ± 0.016) 


0.185 ±0.017 


245 


-1.640 ±0.022 


0.149 ±0.022 


7 


< 0.1 




ACEFGH 


Yes 


377 


(-0.232 ± 0.018) 


0.199 ±0.019 


210 


-1.642 ±0.024 


0.153 ±0.023 


8 


< 0.1 




ACEFG 


Yes 


285 


-0.249 ± 0.021 


0.206 ±0.023 


156 


-1.650 ±0.025 


0.140 ±0.025 


9 


< 0.05 




ACEFGH 


No 


716 


(-0.240 ± 0.014) 


0.186 ±0.014 


226 


-1.651 ±0.024 


0.150 ±0.024 



Table 3. Results of the simulations 



Model 


Ml 


0"RC,I 


A/k 


CRC.K 


Ml (true) 


Mn (true) 


2 


(-0.290 ± 0.010) 


0.240 ±0.007 


-1.648 ±0.018 


0.251 ±0.015 




-1.486± 0.034 


3 


-0.327 ± 0.012 


0.222 ±0.010 


-1.697 ±0.022 


0.238 ± 0.021 


-0.212± 0.020 


-1.485± 0.041 


4 


-0.290 ± 0.023 


0.190 ±0.021 


-1.625 ±0.035 


0.202 ± 0.025 


-0.239± 0.033 


-1.599± 0.043 


5 


-0.264 ± 0.053 


0.214 ±0.040 


(-1.598 ±0.110) 


0.242 ± 0.099 


-0.238± 0.062 




6 


(-0.333 ± 0.015) 


0.218 ±0.011 


-1.702 ±0.025 


0.234 ± 0.021 




-1.538± 0.033 


7 


(-0.337 ± 0.018) 


0.208 ±0.015 


-1.707 ±0.027 


0.233 ± 0.029 




-1.535± 0.036 


8 


-0.330 ± 0.021 


0.217 ±0.018 


-1.715 ±0.030 


0.231 ±0.029 


-0.189± 0.034 


-1.535± 0.039 


9 


(-0.306 ±0.011) 


0.210 ±0.009 


-1.664 ±0.028 


0.228 ± 0.021 




-1.587± 0.037 




^ ° Fig. 3. Fit of Ml versus metallicity and (F — /) colour for model 

Fig. 2. The fit of Eq. [TO]to the data in the J-band (upper panel), 8, with best fits indicated by the dashed line. The slope of the fit 
and K-hwd for model 3 . against colour is not significant. 



lease and increasing the sample to 24 clusters. Van Helshoecht & 

Groenewegen (2007; hereafter vHG) find a value fainter by 0.05 Grocholski & Sarajedini (2002) and vHG used averages over the 
mag, or A/k = — 1 .57±0.05 on the Bessell & Brett system. Both cluster sample to arrive at the quoted means. 
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Fig. 4. Fit of Afx versus metallicity and {V — K) colour for 
model 8, with best fits indicated by the dashed line. The slope of 
the fit against metallicity is not significant. 



To compare these values to the absolute magnitude of the lo- 
cal Hipparcos sample, a "population correction" has to be made, 
to account for the difference in metallicity and age of the RC 
population. The calculation of this correction is outlined and tab- 
ulated in Girardi & Salaris (2001) and Salaris & Girardi (2002). 
In the K-bwd, Salaris & Girardi (2002) show that the correction 
is a strong function of age for ages below 3 Gyr The top panel 
in their Fig. 3 shows that the correction is well behaved for ages 
above about 4 Gyr, and for each metallicity a linear function was 
fitted. 

From vHG the five clusters older that 4 Gyr were taken and 
the population correction was determined by linear interpolation 
in age and [Fe/H] in the results in Salaris & Girardi. The er- 
ror in AMk is based on a 1 Gyr error in age and 0.1 dex in 
metallicity. The results are listed in Table |4] which lists the age 
and metallicity (see vHG for the references for age and metal- 
licity determination), the A/k value and its error that takes into 
account the error in reddening and the assumed distance to the 
cluster (from vHG), the population correction with error, and the 
corrected value. The weighted mean of the corrected A/k values 
is — 1.39± 0.06 with a dispersion of 0.2 mag. Correcting to the 
2MASS system (Cai-penter 2001) this becomes -1.43± 0.06. 
Although the dispersion is large, the mean value is also fainter 
than the straight mean for the cluster sample and agrees within 
2-sigma with the determination from the Hipparcos sample. 

To settle the issue on the importance of the bias and the abso- 
lute X-magnitude of RC stars would require accurate NIR mag- 



nitudes of a 100 to a few hundred (cf. Table 2) bright {K<h) RC 
stars. Given the brightness, this represents a challenge to modern 
instrumentation because of saturation. 

The absolute magnitudes derived here are in better agree- 
ment with theory than previously. Using a plausible star for- 
mation history for the solar neighbourhood, Salaris & Girardi 
(2002) derived absolute magnitudes of Mi — —0.17 and A/k = 
— 1.584 (on the 2MASS system), which they compared to the 
Hipparcos-based results quoted in Alves et al. (2002) of Afj — 
-0.26 ± 0.03 and A/k = -1.644 ± 0.03 (on the 2MASS sys- 
tem), which indicates differences at the 2-3 sigma level. The new 
results of A/j -0.22 ± 0.03 and A/k = -1.54 ± 0.04 agree 
with theory at the 1-2 sigma level. 

Finally, some implications for existing distance determina- 
tions using the RC are discussed. The derived value of the ab- 
solute magnitude in the /-band of A/j — —0.22 ± 0.03 is not 
very different from previously adopted values in the literature 
of —0.23 or —0.26, so the impact on derived distances is not so 
great. 

The distance to the LMC quoted by Pietrzynski & Gieren 
(2002) of 18.501 ± 0.049 (random-nsystematic error and assum- 
ing no population correction) based on JK observations of two 
fields in the bar would become 18.35 ± 0.05 for an assumed 
population correction of —0.03 (Salaris & Girardi 2002). 

Babusiaux & Gilmore (2005) obtained infrared data on some 
fields in the direction of the Galactic bulge. Their distance scale 
is tied to one field for which they obtain A/k = —1.72 for an 
assumed distance of 8 kpc. Taking Myl — —1.54 would result in 
a distance of 7.4 kpc (±0.2 kpc judging from the information in 
their Table 2). Adding the error in A/k results in a total error of 
0.3 kpc. 

Using similar type of observations Nishiyama et al. (2006) 
obtained 7.52 ± 0.1 kpc for A/k = -1.59 ± 0.03, which trans- 
lates to 7.33 ± 0. 14 kpc with the new DM. 

This can be compared to a distance of 7.94 ± 0.37 (random) 
± 0.26 (systematic) kpc from Groenewegen et al. (2007) based 
on /C-band monitoring observations of Population-Il cepheids 
and RR Lyrae stars (distance tied to an assumed LMC distance 
modulus (DM) of 18.50) and the direct determination of 7.73 
± 0.32 kpc based on the orbit of the star called S2 around the 
central black hole when post-Newtonian physics is taken into 
account (Eisenhauer et al. 2005, Zucker et al. 2006), although 
neglect of the space motion of Sgr A* may result in systematic 
errors of about 0.1-0.45 kpc (Nikiforov 2008). All these inde- 
pendent observation are consistent and seem to indicate a DM to 
the GC of slightly less than 8 kpc. 
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